Disorder-Induced First Order Transition and Curie Temperature Lowering in 

Ferromagnatic Manganites. 
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We study the effect that size disorder in the cations surrounding manganese ions has on the 
magnetic properties of manganites. This disorder is mimic with a proper distribution of spatially 
disordered Manganese energies. Both, the Curie temperature and the order of the transition are 
strongly affected by disorder. For moderate disorder the Curie temperature decreases linearly with 
the the variance of the distribution of the manganese site energies, and for a disorder comparable to 
that present in real materials the transition becomes first order. Our results provide a theoretical 
framework to understand disorder effects on the magnetic behavior of manganites. 

PACS numbers: 75.47.Gk,75.10.-b. 75.30Kz, 75.50.Ee. 
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Compounds with chemical formula R^^A^MnOa 
(R=trivalent element and A= divalent element) are usu- 
ally called manganites. These materials exhibit an ex- 
traordinarily large (so called colossal) magnetoresistance. 
The study of the origin of the colossal magnetoresis- 
tance as well as its possible technological applications 
has triggered in the last years an intensive research on 
these materials pj, |3, This research has revealed a rich 
variety of exotic phenomena, and after twenty years of 
intense study new and exciting phenomena still appear 
when studying manganites |4|, |5|, |6|, 13 • 

In manganites the ferromagnetic order is driven by 
the motion of the carriers and, therefore, their proper- 
ties depend on the competition between kinetic energy, 
tending to delocalize the carriers, and localization ef- 
fects, such as antiferromagnetic coupling between the 
Mn core spins and the Jahn- Teller coupling. Disorder 
also reduces the carriers mobility, and, thus, it strongly 
affects the stability of the ferromagnetic (FM) phases. 
Recently it has been possible to synthesize cation or- 
dered Lno.sBao.sMnOa, where Ln is a rare earth and LnO 
and BaO planes alternates along the c axisQ- These 
half-doped manganites may give a good benchmark to 
test theories of the influence of disorder on the elec- 
tronic and magnetic properties of manganites. Experi- 
mental work in FM Lao.sBao.sMnOa shows[j| that dis- 
order strongly suppress the transition temperature and 
changes the character of the transition, form continuous 
to weakly first order. 

Attficld and coowor kers[l3[lll found that in FM per- 
ovskites of the form AMnC>3, the variation of Curie tem- 
perature with disorder is related with the distribution of 
the A cations radius, r^, and there is not significant de- 
pendence on the A site charge variance. Therefore, we 
expect that the observed 9] difference in T c between the 
ordered and disordered Lao.sBan.sMnOa perovskites is 
related with the different distributions of ta- Note that 
standard ionic radii values are ta — 1-21 for La 3+ and 
r A = 1.47 for Ba 2+ llj. In the ordered Lao. 5 Bao.5Mn0 3 
perovskite, Fig^a), alternating planes of LaO and BaO 
separate Mn02 sheets in such a way that each Mn ion 
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FIG. 1: (color online) Schematic representation of an ordered 
(a) and disordered (b) Lao.sBao.sMnOa perovskite. In the 
ordered case planes containing La and Ba cations alternate 
along the z-direction. In the disorder case the cations are 
randomly distributed. 



has four La 3+ neighbors in a direction and four Ba 2+ 
neighbors in the opposite direction being all the Mn ions 
equivalent. In the disorder perovskites, Fig^b), the La 
and Ba cations are randomly distributed and the Mn's 
are surrounded by different combinations of divalent and 
trivalent cations. These different chemical environments 
traduce in different chemical pressures on the Mn ions, 
and is a source of disorder in the Mn electronic levels. 

In this work we combine Ginzburg-Landau formalism 
with realistic two orbital DE microscopic calculations, 
to study the effect that the disorder in the position of 
the trivalent and divalent cations has on the magnetic 
properties of FM manganites. The main results of the 
work are, i) The Curie temperature, T c decreases with 
the disorder, ii) as the disorder increases, the FM-PM 
phase transition becomes more abrupt and eventually, 
for a critical disorder, the magnetic transition becomes 
first order and iii) we identify the variance of disorder 
distribution as the relevant parameter for the decrease of 
T c . All these three points agree with the experimental 
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results H IHIT3- 

Our analytical results provide a frame- 
work for understanding the dependence of T c on disorder, 
both, in experiments and in Monte Carlo simulations of 
simplified models HI Q E 13 . 

Model. In a manganite of formula A^MnOs, there 
are A — x electrons per Mn ion. The crystal field splits the 
Mn d levels into an occupied strongly localized ti g triplet 
an a doublet of e g symmetry. Furthermore, the Hund's 
coupling in Mn ions is very large and aligns the spins of 
the d orbitals. Effectively, there are 1 — x electrons per 
Mn ions hoping between the empty e g Mn states. 

The large Hund's coupling forces each electron spin 
to align locally with the core spin texture. The spin 
of the carriers is conserved in the hoping between Mn 
ions, being the tunnelling amplitude maximum and the 
kinetic energy minimum, when the spins of the Mn ions 
are parallel and the system is ferromagnetic. This is ba- 
sically the so called double exch ang e (DE) mechanism 
proposed fifty years ago by ZenerJlTj to explain the FM 
order in manganites. When the temperature increases 
the kinetic energy minimized by the FM order competes 
with the orientational entropy of the Mn core spin and 
at T c , the system becomes paramagnetic (PM)|lg- I n 
the PM phase the reduction of the kinetic energy could 
favor that localization effects become more effective and 
a metal insulator transition could occur near the FM-PM 
transition[19j. In the case of Lao.sBao.sMnOs, the local- 
ization effects are rather weak and the material is metal- 
lic at both sizes of the FM-PM transition. Therefore, in 
our model we do not consider coupling with the lattice, 
however as the DE mechanism depends strongly in the 
kinetic energy gain, we treat the motion of the carriers 
in a realistic way by including the two e g orbitals. 

Hamiltonian. With the above considerations, the two 
orbital DE Hamiltonian takes the following form, 

h= £ /,., • (i) 

<z,j>,a,b i, a 

here Cf a creates an electron in the Mn ion located at site 
i in the e g orbital a (<z=l,2 l=|a; 2 — y 2 > and 2=|3z 2 — 
r 2 >). In the limit of infinite Hund's coupling, the spin 
of the carrier should be parallel to the Mn core spin Sj, 
and the tunnelling amplitude is modulated by the spin 
reduction factor 

/. = cos — cos — ^ + e«<*-**) sin ^ sin *L (2) 
Jhj 2 2 2 2 ' 

where cf>i} are the Euler angles of the Mn core spins. 
The hopping amplitude depends both on the direction u 

between sites i and j and the orbitals involved; t^f ■* = 
±Vlf l [l ) = ±VSt^ ] = 3fffi = t. In the z direction the 
only nonzero term is if 2 = [|. Hereafter t is taken 
as the energy unit. 

The last term describes the diagonal disorder. hi >a = 
Cf a Ci^ a is the occupation operator of orbital a at site i 
and 6i is the energy shift produced by the chemical pres- 
sure on the Mn ion at site i. This shift affects equally 



both e g orbitals. As discussed above, the chemical pres- 
sure at site i depends on the ionic radii of the cations 
surrounding the Mn ion. We assume that the total shift 
at site i is the sum of the shifts produced by the eight 
next neighbors cations. Cations with big ionic radius in- 
duce a positive energy shift A whereas cations with small 
ionic radius induce a negative energy shift — A. With this 
election a Mn ion surrounded by equal number of small 
and large cations has zero energy shift. In the case of 
a manganite of composition R 5 A .5MnO3 the diagonal 
shift takes values in the range — 8 A < < 8 A. In this 
description of the cationic disorder the strength of the 
disorder is defined by the value of A. 

Free energy. In order to describe thermal effects we 
have to compute the free energy. In manganites the 
Fermi energy of the carriers is much greater than the 
typical Curie temperature and we only consider the en- 
tropy of the classical Mn core spins. In the mean field 
approximation 20J and in the limit of small magnetiza- 
tion, the entropy per Mn core spin takes the form, 

log2 3 2 9 4 
S(m) = — --m -y Q m + ... (3) 

where m is the thermal average of the relative magneti- 
zation of the Mn core spins. 

At finite temperatures the magnetization is not satu- 
rated and the spin reduction factor, Eq. @ is smaller 
than unity. Treating the spin fluctuations in the vir- 
tual crystal approximation, the spin reduction factor is 
substituted by its expectation value, that for small mag- 
netization has the form PH. 

h^f m = \-\m>-^ + ... (4) 

With this the thermal average of the internal energy per 
Mn ion can be written in the form, 

E (E t u a , b < C+ a C^ b > +± £ e i < h i>a >) (5) 

being N the number of Mn ions in the system. In Eq. 
, it is evident that the relative importance of the disor- 
der increases when the magnetization, and thus the spin 
reduction factor decreases. For small values of A the in- 
ternal energy can be expanded in powers of the disorder 
strength. As the mean value of e, is zero, the first term 
different from zero is proportional to the variance, <r 2 (ei), 
of the diagonal disorder distribution, 

E = f m (^E + a^p^j , (6) 

where Eq < is the kinetic energy per Mn ions in the 
disorder free FM phase, and 

i 

In Eq. ©, the coefficient a is negative as the electrons 
prefer to place on sites with negative values of e^. 
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Combining Eq. il-SITI) . the dependence of the free en- 
ergy, F = E — TS, on the magnetization, can be written 
as 
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£ -yaa 2 (e 4 ) ) n,< + ... (M 



and the Curie temperature takes the form 
4 /„ 9 



Tr 



ibK E a --aa\e t ) 



(9) 



where is clear that T c decreases when the disorder 
strength increases. The order of the transition can be 
inferred from the sign of the quartic term in Eq. (fHJ . If 
the quartic term is positive the transition is second order, 
while a negative quartic term implies the existence of a 
first order phase transition. Using the previous expres- 
sion for T c , a first order transition takes place if, 
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Equations @ and (|10(1 are the main results of this work, 
they show that the Curie temperature decreases when 
the strength of the disorder increases and that for disor- 
der strong enough the FM-PM transition changes from 
second to first order. Comparing Eq. © and Ea. Hl(J|) it 
results that for changes in the Curie temperature larger 
than 10 per cent, the FM-PM transition transforms from 
second to first order. 

Numerical results. In order to check the approxima- 
tions we have done and the validity of the obtained re- 
sults, we diagonalize numerically the ferromagnetic two 
orbital DE Hamiltonian, in presence of diagonal disorder, 
Eq. QJ. We consider the case of disorder in the ionic radii 
of the cations surrounding the Mn ions. As discussed 
above, this disorder produces a distribution of diagonal 
energy shifts in the range — 8A < e, < 8A. For a random 
distribution of the two type of cations, the variance of the 
diagonal disorder distribution is <J 2 {ei) — 8A 2 . We have 
diagonalized the Hamiltonian for different disorder real- 
ization and different values of the disorder strength A. 
We work with a cluster containing 12x12x12 atoms and 
with periodic boundary conditions. By studying smaller 
clusters we have checked that our results are free of finite 
size effects. 

In Fig|2Ji we plot, for the full polarized FM phase, the 
internal energy per Mn ion as function of the disorder 
strength A. As expected 22], for small values of the dis- 
order the internal energy decreases quadratically with A. 
Taking into account the dependence of the internal en- 
ergy on the spin reduction factor, in mean field approxi- 
mation the Curie temperature gets the form 
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FIG. 2: Polynomial fit to a) internal energy and b) Curie 
temperature as obtained from numerical simulations. Bars in 
a) indicates the typical indetermination in the numerical cal- 
culations. In b) the dashed line corresponds to the quadratic 
approximation to T c (Eq. ©). In b) the regions with first 
and second order FM-PM transition are shown. 



In Fig^, we plot T c , as obtained numerically from the 
internal energy, as function of A. In the disorder free 
case, T c ~ 0.8i that compared with the Curie tempera- 
ture of the ordered Lao.5Ba .5Mn0 3 , T C =350K, 0, im- 
plies a value of the hopping amplitude of t ~ 0.038eV. 
For changes up to 20 per cent, we find that the de- 
crease of T c with A can be fitted very accuracy with 
a quadratic dependence. Experimentally the differences 
in T c between the ordered and disordered samples in 
Lao.5Bao.5MnC"3 is near 14 per centQ, therefore we ex- 
pect that for manganites with atomic radii size disorder, 
the expressions obtained analytically for T c , Eq. JJjJ, and 
the order of the transition, Eq. (|10fl should be valid. 

In figEb, we also indicate the disorder strength for 
which the magnetic transition changes form second to 
first order. As discussed above, the transition becomes 
first order for changes in T c larger than 10 per cent. 
In the case of Lao.5Bao.5Mn O3, the disorder produces 
changes in T c of near 14 per cent and therefore our re- 
sults explain the experimental observed change in the or- 
der of the magnetic transition^. We expect that in FM 
manganites with smaller differences in the atomic radii 
than in Lao.5Bao.5Mn03, for example Lao.sSrn^MnOs 
(r A = 1.21 for La 3+ and r A = 1.31 for Sr 2+ |l|l) the 
effect of the disorder in T c should be smaller and the 
order of the magnetic transition should be second order 
independently of the order in the position of the cations. 

Finally we want to check that, as proposed by 
Rodriguez-Martinez and Attfieldjll]|, the relevant figure 
for quantifying the disorder is the variance of the disorder 
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FIG. 3: (color online) Curie temperature versus variance of 
diagonal disorder for different models of disorder. The differ- 
ent kinds of disorder are explained in the text 



distribution. We have repeated the numerical calculation 
of T c for different distribution of diagonal disorder. Apart 
from the cation disorder distribution , we have also ana- 
lyzed a uniform distribution of the diagonal energy shift 
between two values, a binary distribution and a gaussian 
distribution of the diagonal disorder. 

The results are plottted in Fig|3 For all models stud- 
ied, we find that for moderate disorder strength, the 
Curie temperature decreases linearly with the variance 
of the disorder. There are deviations from this depen- 



dence in the limit of strong disorder. In this limit the 
carriers start to localize in the sites with smaller energy 
and perturbation theory becomes not valid. This regime 
is not the relevant one in FM manganites, and the re- 
sults shown in Fig. [3] agree with the experimental results 
obtained by Rodriguez-Martinez and AttfieldpT.]. The 
numerical results also support our analytic findings that 
for moderate disorder strength the variation of the Curie 
temperature and the order of the transition depend on 
the the variance of the distribution of the diagonal energy 
shifts in the electronic Hamiltonian. 

In summary, we have studied a realistic model of man- 
ganites by means of exact microscopic calculations and 
Landau Theory formalism. The different sizes of cations 
surrounding a manganese ion are included in the model as 
an energy shift in that manganese site. This model repro- 
duces the experimental results. It explains the observed 
strong reduction of T c in disordered samples with respect 
to ordered ones[2j. Moreover, it demonstrates, in agree- 
ment with experiments 0, that disorder makes the FM- 
PM transition more abrupt and, for a enough disorder 
strength, this transition becomes first order. Finally our 
formalism identifies in a natural way the variance of the 
distribution of the diagonal energy shifts, as the relevant 
parameter to characterize changes in the ma gnetic p rop- 
erties, independently of the model of disorder [ill Il2l . l28| . 
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